Draft version February 7, 2013 

Preprint typeset using L^T^X style emulateapj v. 5/2/11 



m 
o 



Oh! 

o 
is 

c3 



HEATING OF THE MAGNETIC CHROMOSPHERE: OBSERVATIONAL CONSTRAINTS FROM CA II A8542 

SPECTRA 

J. de la Cruz Rodriguez 1 , B. De Pontieu 2 , M. Carlsson 3 , and L.H.M. Rouppe van der Voort 3 

1 Department of Physics and Astronomy, Uppsala University, Box 516, SE- 75120 Uppsala, Sweden 
2 Lockheed Martin Solar & Astrophysics Lab, Org. A021S, Bldg. 252, 3251 Hanover Street Palo Alto, CA 94304 USA and 
institute of Theoretical Astrophysics, University of Oslo, P.O. Box 1029 Blindern, N-0315 Oslo, Norway 

Draft version February 7, 2013 

ABSTRACT 

The heating of the Sun's chromosphere remains poorly understood. While progress has been made 
on understanding what drives the quiet Sun internetwork chromosphere, chromospheric heating in 
strong magnetic field regions continues to present a difficult challenge, mostly because of a lack of ob- 
servational constraints. We use high-resolution spectropolarimetric data from the Swedish 1-m Solar 
Telescope to identify the location and spatio-temporal properties of heating in the magnetic chromo- 
sphere. In particular, we report the existence of raised-core spectral line profiles in the Ca II A8542 
line. These profiles are characterized by the absence of an absorption line core, showing a quasi-flat 
profile between A « ±0.5 A, and are abundant close to magnetic bright-points and plage. Comparison 
with 3D MHD simulations indicates that such profiles occur when the line-of-sight goes through an 
"elevated temperature canopy" associated with the expansion with height of the magnetic field of flux 
concentrations. This temperature canopy in the simulations is caused by ohmic dissipation where 
there are strong magnetic field gradients. The raised-core profiles are thus indicators of locations of 
increased chromospheric heating. We characterize the location and temporal and spatial properties 
of such profiles in our observations, thus providing much stricter constraints on theoretical models of 
chromospheric heating mechanisms than before. 

Subject headings: Sun: chromosphere — Sun: magnetic topology — Sun: faculae, plages — line: 
formation — magnetohydrodynamics — polarization 
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1. INTRODUCTION 

The heating of the Sun's chromosphere requires 
more than an order of magnitude more mechan- 
ical energy flux than the corona and heliosphere 
combined (Anderson & Athay 1989; Withbroe & Noyes 
1977). Nevertheless, the heating mechanism powering 
the chromosphere remains elusive, especially in the mag- 
netic chromosphere in and around network and plage 
regions (e.g., Judge et al. 2010). This is in part be- 
cause most studies have focused on quiet- Sun internet- 
work regions (e.g., Carlsson & Stein 1997; Straus et al. 
2008). In addition, most theoretical studies have avoided 
detailed comparisons with chromospheric diagnostics, 
which are most often formed under non-LTE and non- 
equilibrium conditions and thus difficult to model, in- 
stead focusing on simply comparing the energy flux as- 
sociated with specific physical heating mechanisms (e.g. 
Goodman & Kazeminezhad 2010) to the canonical val- 
ues derived from semi-empirical ID hydrostatic models of 
the chromosphere (Vernazza et al. 1981; Fontenla et al. 
1993). Such a limited comparison does not properly 
capture the chromospheric conditions, which are sig- 
nificantly different from the ID VAL and FAL mod- 
els (in terms of dynamics and structuring, see, e.g., 
Carlsson & Stein 1997). 

Here we report on new observational constraints on 
chromospheric heating derived from the Ca II A8542 
spectral line. During the past decade, the Ca II in- 
frared (IR) triplet lines (A8498, A8542, A8662) have be- 
come popular chromospheric diagnostics in solar stud- 
ies. In particular, the Ca II A8542 line has been 
used in a variety of different chromospheric stud- 



ies, for example, to determine the disk counterpart 
of spicules (Rouppe van der Voort et al. 2009), analyze 
the energy flux carried by acoustic waves into the 
lower atmosphere (Vecchio et al. 2009), study torsional 
motions in the atmosphere (De Pontieu et al. 2012; 
Wedemeyer-Bohm et al. 2012), or to carry out spec- 
tropolarimetric studies in sunspots (Socas-Navarro et al. 
2000; Lopez Ariste et al. 2001; Kleint 2012). 

Cauzzi et al. (2009) suggested a connection between 
chromospheric heating and the width of the H I A6563 
(Ha) line, which is found to be strongly sensitive to 
changes in temperature. Furthermore, their maps of the 
line core intensity of Ca II A8542 show strong similarities 
with their Ha line- width maps. Further support for 
a relationship between core brightness of Ca lines and 
elevated temperatures is provided by Henriques (2012), 
who used the intensity of Ca II A3968 as a measurement 
of the upper-photospheric temperature. 

Although the velocity field and temperature stratifica- 
tion in the solar atmosphere produce an imprint in the 
spatially-resolved profiles from the Ca II IR lines, they 
usually remain absorption lines with a conspicuous line 
core. Interestingly, Pietarila et al. (2007b) report some 
examples of reversals in the NLTE core of the Ca II A8542 
line, which are later associated with elevated chromo- 
spheric temperatures in Pietarila et al. (2007a). However 
the mechanism producing this increase in temperature 
or the spatio-temporal properties are not investigated. 
These raised-core (RC hereafter) profiles are found in 
pixels with intermediate magnetic flux in their observa- 
tions. 

In this letter, we use full-Stokes Ca II A8542 data to 
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Figure 1. Ca II A8542 filtergrams (top-row): Stokes I at —1055 mA from line center (left), Stokes V at the same wavelength (middle) 
and the line core intensity (right). The color markers indicate the location of hand-picked quiet-Sun profiles (green), RC profiles (orange) 
and bright-point profiles (blue). The circular markers indicate the location of the intensity profiles shown in each of the panels below (same 
color coding). 



show the ubiquity of such RC profiles in the vicinity of 
strong magnetic fields. We report on the peculiar spatial 
distribution of the RC profiles, and analyze their spatial 
and temporal properties to provide constraints on chro- 
mospheric heating models. Finally, we propose a scenario 
that explains the profile shapes, using synthetic spectral 
line profiles from a 3D radiative MHD simulation. 



2. OBSERVATIONS 

The observations are of a bipolar active region 
(AR10998) at approximately 205" heliocentric distance. 
The datasets were acquired at the Swedish 1-m Solar 
Telescope (SST) using the CRisp Imaging SpectroPo- 
larimeter (CRISP, Scharmer et al. 2008) on 14 Jun 2008 
at 08:36:54. 

Our data consists of a 41 minute time-series including 
225 full-Stokes line-scans (time steps) of the Ca II A8542 
line, with a temporal cadence of 11 seconds. The line was 
sampled at 29 line positions across the range ±1.94 A 
from the line center, with a spectral sampling of 97 mA 
in the core up to ±0.78 A, and 194 mA in the outer 
wings. 

The data are processed using the image restoration 
technique Multi- Ob ject-Multi- Frame-Blind- Deconvolu- 
tion (MOMFBD, van Noort et al. 2005), which allows to 
achieve a spatial resolution of X/D =0 / /18. Instrumen- 
tal polarization is compensated using the 2D calibration 
described in Schnerr et al. (2011) and the SST model at 
854 nm from de la Cruz Rodriguez (2010). 

We have made extensive use of CRISPEX 
(Vissers & Rouppe van der Voort 2012) a tool for 
the exploration of CRISP datasets. 



2.1. The raised- core profiles 

The observed field-of-view (FOV) contains two patches 
of network with opposite polarities. In the photosphere, 
observed in the wings of the A8542 line, the images show 
granulation with two concentrations of bright-points. In 
the chromosphere, visible in the core of the line, the im- 
ages show a landscape of dark elongated fibrils that ex- 
tend outwards from the network patches. 

A section of our FOV is enlarged for clarity in Fig. 1. 
This patch of mostly unipolar plage shows strong mag- 
netic field concentrations (Stokes V top middle panel) 
that are tightly correlated with the locations of bright- 
points (Stokes /, top left panel). The line core image (top 
right panel) , which is sensitive to chromospheric plasma, 
shows fibrils extending from the plage region. The plage 
region is generally very bright compared with the fibrils. 

We distinguish three different "typical" line profiles and 
give examples from three different locations for each pro- 
file type: quiet (locations and spectra marked with green 
symbols), bright-point (blue symbols) and RC profiles 
(orange symbols). The quiet profiles show a "classical" 
Ca II A8542 profile with bright photospheric wings and a 
chromospheric absorption core (green profiles in bottom 
row of Fig. 1). The shape of the bright-point profiles 
(blue in bottom row) is very similar to that of the quiet- 
Sun profiles but with enhanced intensity for all wave- 
lengths. 

The RC profiles look markedly different. The pho- 
tospheric wings of RC profiles are virtually identical in 
brightness and shape to those of quiet profiles. However, 
the chromospheric core of RC profiles is very different: 
it is almost flat, and much brighter than the core of the 
quiet profiles. The chromospheric core of RC profiles ac- 
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Figure 2. Spatial distribution of RC profiles measured as number of zero-crossings of d//dA over the core region (±582 mA). The left 
panel shows a Stokes V map integrated over all 29 line positions with contours delineating areas with 3 (red) and 4 (blue) zero-crossings. 
The rightmost panels show a zoom-in on the region outlined in the left panel, with a wide-band image as background at the top. An online 
movie is available. 



tually has an intensity that is similar to that of the core 
of the bright-point profiles. 

In a large fraction of the RC profiles, we have also 
noticed two weak emission lobes at approximately A = 
±194 mA from line center. These lobes are further dis- 
cussed in §§ 2.3 and 3. 

2.2. Spatial distribution and ubiquity 

One clue to the cause of the RC profiles comes from 
investigating their spatial distribution and ubiquity in 
the plage region. Clearly not all locations show RC pro- 
files. To determine, in an automated fashion, the loca- 
tion of RC profiles we tried several approaches. We found 
that the ratio between core and wing intensities reveals 
many locations of raised core profiles, but visual inspec- 
tion shows that it also flags many "false positive" loca- 
tions (usually narrow absorbing cores that are strongly 
Doppler-shifted). A better method exploits the fact that 
around the core of these profiles the intensity is relatively 
flat with wavelength, so that the derivative of intensity 
with wavelength d//dA hovers around 0, with frequent 
changes of the sign of the derivative. We thus developed 
an algorithm that calculates the number of zero-crossings 
of d//dA around the core of the line. Both bright-point 
and quiet profiles have values of 1, whereas RC profiles 
almost always have larger values. 

By overplotting contours of the number of zero- 
crossings of d//dA (red and blue contours in Fig. 2), we 
find that the RC profiles are located in the proximity 
of photospheric bright-points and in extended patches of 
photospheric granulation, which are embedded in patches 
of plage. In fact, they are quite rare outside this config- 
uration. This is quite clear from the enlarged regions 
shown in the right column of Fig. 2. The RC profiles 



are preferentially located inside the ring of bright-points, 
with only a few exceptions found outside. Our analy- 
sis also shows that these RC profiles are ubiquitous and 
occur wherever there is plage (or network), with their 
frequency increasing in more complex regions. The spa- 
tial extent of seemingly coherent or contiguous regions 
with RC profiles is of order a fraction of an arcsecond up 
to an arcsecond. 

2.3. Temporal properties 

The temporal behavior of the RC profiles is very dif- 
ferent from that of a more typical profile, as can be 
seen by comparing the three top-right panels with the 
second panel from the left in the top row of Fig. 3. 
The latter shows the temporal evolution of the Ca II 
A8542 line profile in a location that does not show 
strong RC behavior. It is dominated by a succes- 
sion of magneto- acoustic shocks, visible as strongly 
absorbing features and similar to those described by 
Langangen et al. (2008); Vecchio et al. (2009). Two brief 
moments of RC profiles occur in between shocks. The 
three RC locations (rightmost panels, top row) on the 
other hand show raised core emission for much of the 
time-series, with only occasional shock-induced absorp- 
tion. The RC profiles are not constant in time but show 
significant variability. We find that the imprint from 
the red emission lobe at ±194 mA is often noticeable 
throughout the timeseries, which is not always the case 
for the blue lobe. 

The online movie of Fig. 2 shows evidence of contigu- 
ous and/or coherent spatial regions that contain RC pro- 
files with typical spatial and temporal scales of order 
0.5"and 1-2 minutes. These regions change shape and 
move about within the plage region. To determine the 
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Figure 3. Temporal evolution and lifetime of the RC profiles. The top row illustrates the temporal evolution of some A8542 profiles at 
locations indicated with color markers in the top left panel and repeated in each A — t plot (right panels). The second panel from left 
corresponds to a traditional absorption profile showing a shock induced temporal pattern. The last three panels correspond to the time 
evolution of three RC profiles. The time of the image in the top-left panel is indicated in the time slices with a yellow marker. The bottom 
row provides information about the typical timescales that dominate the temporal evolution of the RC profiles: map of the auto-covariance 
timescale (8x rebinned) (left) and a histogram of auto-covariance timescales for locations with RC profiles (right), using the original data 
(dashed-red) , 8x8 (solid-black) and 16 x 16 rebinning (dotted-blue) . The vertical black-dotted line indicates the median auto-covariance 
timescale: 81 seconds. 



typical timescales over which RC profiles change, we com- 
puted, for each RC location, the auto-covariance func- 
tion of a timeseries that is equal to 1 when the number 
of zero-crossings is larger than 1, and otherwise 0. The 
auto-covariance timescale for that location is then the 
full width half max of this function. 

To make sure our timescale reflects the intrinsic life- 
time and is not dominated by either seeing-induced or 
proper motion of the coherent RC region within the FOV, 
we calculated the auto-covariance timescale for data that 
was rebinned with a range of different parameters (from 



2x2tol6xl6 pixels). The timescales increase as the 
rebinning parameters are increased from lxltol6xl6 
(Fig. 3, lower right panel). This is because the timescales 
of the data at the original resolution are dominated by 
seeing deformations, jitter and proper motions leading to 
short lifetimes. Once we rebin the data 8x8 (solid-black) 
we no longer see a difference in timescale compared to 
16 x 16 (dotted-blue) rebinned data, with both showing 
median values of about 80 seconds. This means that the 
intrinsic dominant lifetime of these coherent RC regions 
is about 1.5 minutes. It also implies that the typical 
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Figure 4. Synthetic A8542 observations and physical quantities from our 3D MHD simulation. In the top row, monochromatic intensity 
in the wing (left), the vertical component of the magnetic field at z = km (middle) and the monochromatic intensity at line center (right). 
The lower-left panel shows the intensity profiles at locations indicated with colored markers in the top row: quiet-Sun (green), bright-point 
(blue) and raised-core profile (orange). The lower-middle panel shows the temperature stratification (dashed) and source function (solid) 
at each location. The markers over-plotted on the source function indicate the height at which r v = 1 in the continuum (plus) and at line 
center (star). The lower right panel shows a 2D cut (along the orange line shown in the top middle panel) of the temperature as a function 
of height in the simulation. The zero-point in x is at the location of the orange circle in the top middle panel, the site of an RC profile. 
Markers show the location of r v = 1 in the continuum (plus) and at line center (star). 



spatial scales are less than 8 pixels, i.e., 8 x 0.07 ~ 0.5". 

3. 3D MHD SIMULATION 

We use a snapshot from a 3D MHD simulation cal- 
culated with the Bifrost code (Gudiksen et al. 2011). 
The simulation has an extent of 24 x 24 x 16.8 Mm, 
encompassing the upper convection zone, photosphere, 
chromosphere and corona. In the vertical direction, 
the simulation extends from 2.4 Mm below to 14.4 Mm 
above average optical depth unity at A = 500 nm. The 
grid spacing is 48 km horizontally and non-equidistant 
vertically with a spacing of 19 km between z = — 1 
and 5 Mm and increasing towards the upper and lower 
boundaries. The simulation includes optically thick ra- 
diative transfer including scattering in the photosphere 
and low chromosphere, parametrized non-LTE radiative 
losses in the upper chromosphere, transition region and 
corona, thermal conduction along magnetic field lines 
and an EOS that includes the effects of non-equilibrium 
ionization of hydrogen, see Gudiksen et al. (2011) and 
Carlsson & Leenaarts (2012) for details. The magnetic 
field in the simulation consists of two patches of different 
polarity separated by 8 Mm (Fig. 4, top middle panel). 
The field is rather weak compared with the plage-like re- 
gions of the observations; the average unsigned flux in 
the photosphere is 50 G. The same snapshot was used in 
Leenaarts et al. (2012). 

The synthetic A8542 profiles are computed using the 
code Multi (Carlsson 1986), assuming each vertical col- 
umn in the 3D snapshot to be an independent ID plane- 



parallel atmosphere. The profiles have been convolved 
with the CRISP spectral PSF (FWHM of 111 mA). 

Our synthetic profiles show evidence of quiet, bright- 
point and RC profiles (Fig. 4) that are very similar to 
the observations. By comparing the physical variables 
and radiative transfer properties for various locations and 
profiles, we can investigate the physical cause for the 
shape of the latter. 

At each wavelength, the intensity is approximately 
equal to the source function at monochromatic optical 
depth unity (the Eddington-Barbier relation). The lower 
middle panel of Fig 4 shows the source function at the 
three chosen locations with optical depth unity marked 
for the continuum and for the line center. The inten- 
sity in the wings of the bright-point profile is high be- 
cause the strong magnetic flux regions are evacuated 
compared to their surroundings to stay in horizontal to- 
tal pressure balance. As a result, we find lower densities 
and therefore a deeper formation (Steiner et al. 2001; 
Schiissler et al. 2003; Shelyag et al. 2004; Carlsson et al. 
2004) where the temperature is higher (plus-signs in the 
lower- middle panel). Both the bright-point and quiet- 
Sun locations have monotonically decreasing source func- 
tions with height, giving pure absorption lines, because 
the source function decouples from the Planck function 
lower in the atmosphere than the temperature increase. 
At the RC location we go from a quiet-Sun temperature 
to an increased temperature typical of a magnetic region 
when we cross the expanding magnetic field at about 
0.5 Mm height (see also lower-right panel). The den- 
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sity is still high enough to ensure a coupling between the 
Planck function and the source function and we get a lo- 
cal maximum of the latter, leading to two weak emission 
peaks in the RC profile at roughly ±200 mA. Therefore, 
non-LTE effects and chromospheric heating can explain 
the presence of those emission peaks. The line core is 
formed above the local source function maximum and 
there is therefore a central reversal. The optical depth 
unity at line center is at equal heights for the bright- 
point location and the RC location (because both lo- 
cations are basically magnetic above the canopy height) 
and the source functions are also almost the same, giving 
similar line core intensity at the two locations. 

What makes RC profiles special is that they reveal the 
locations of a sharp temperature gradient with height 
between a relatively quiet photosphere and high temper- 
ature chromosphere. In the simulations this occurs in the 
immediate vicinity of bright-points where we see the hot 
canopy of magnetic field that expands with height from 
its roots in the bright-points (lower-right panel, Fig. 4). 

Our simulations also help explain why our detection 
method of zero-crossings works well in identifying RC 
profiles, given that the two emission lobes are naturally 
produced by a hot canopy overlying the quiet photo- 
sphere. 

4. DISCUSSION AND CONCLUSIONS 

In this letter we characterize the observational prop- 
erties of profiles of the Ca II A8542 line in which the 
core emission is raised compared to typical, quiet, pro- 
files. These RC profiles have photospheric wings similar 
to those of quiet profiles, but line-core intensities that 
are as strong as in the profiles from photospheric bright- 
points. The RC profiles are ubiquitous throughout plage 
and network regions and predominantly located in the 
surroundings of bright-points, usually on the inside of 
the plage region. We find that regions of RC profiles 
show coherent spatial and temporal behavior on scales 
of O^'and 1.5 minutes, respectively. 

We use RC profiles in synthetic observations computed 
from a 3D radiative MHD simulation to determine what 
causes the peculiar line profile. The simulation shows 
that RC profiles occur for locations with a steep increase 
in the temperature stratification of about 1500 K, be- 
tween z ~ 0.5 Mm and z ~ 1.0 Mm. The reduced line 
opacity (due to the vacuum effect of magnetic pressure) 
has a contributive effect to the line core brightening. 

In the simulations the strong temperature gradient 
occurs because the line of sight jumps from a quiet 
photosphere dominated by granulation to the hot chro- 
mospheric canopy associated with flux concentrations. 
These concentrations expand with height from their pho- 
tospheric bright-point roots to a volume filling field 
higher up. Chromospheric heating is caused by strong 
field gradients and currents (and associated dissipation) 
in the vicinity of flux concentrations. This scenario is 
fully compatible with our observations that show RC pro- 
files in the immediate vicinity (but not on top) of bright- 
points, and can explain the enhanced line core intensity 
over the entire plage region. 

Our observations provide strict constraints on the lo- 
cation and spatio-temporal properties of chromospheric 
heating in magnetic regions. The observations fit well 
with the chromospheric heating present in our simula- 



tions, which is dominated by current dissipation. While 
the spatial scale of the dissipation is mostly set by nu- 
merical resistivity in the simulation, recent work suggests 
that ion-neutral interactions lead to a Pedersen resistiv- 
ity that has the same order of magnitude as the numerical 
resistivity, thus rendering these simulations surprisingly 
close to solar conditions ( Martinez- Sykor a et al. 2012). 
Our observations of apparently strong heating at the in- 
terface of interacting flux tubes are also compatible with 
the current-driven heating in the Bifrost models. Fu- 
ture work will have to clarify which of the several other 
physical mechanisms proposed to drive heating in the 
magnetic chromosphere is compatible with the observed 
properties, whether it is high frequency wave heating 
(Hasan & van Ballegooijen 2008; Vigeesh et al. 2009), 
reconnection related to weak granular fields (Lites et al. 
2008; Isobe et al. 2008), or a turbulent cascade of Alfven 
waves (van Ballegooijen et al. 2011). 

Finally, we note that the rather flat line core of the RC 
profiles means that the weak-field approximation (where 
the field strength is inversely proportional to d//dA) 
breaks down, leading to strong noise in the derived mag- 
netic field values. 
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